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ABSTRACT 

We present millimetre (ram) and submillimetre (submm) photometry of a sample of 
five host galaxies of Gamma Ray Bursts (GRBs), obtained using the MAMB02 and 
SCUBA bolometer arrays respectively. These observations were obtained as part of 
an ongoing project to investigate the status of GRBs as indicators of star formation. 
Our targets include two of the most unusual GRB host galaxies, selected as likely 
candidate submm galaxies: the extremely red (i? — if « 5) host of GRB 030115, 
and the extremely faint {R > 29.5) host of GRB 020124. Neither of these galaxies is 
detected, but the deep upper limits for GRB 030115 impose constraints on its spectral 
energy distribution, requiring a warmer dust temperature than is commonly adopted 
for submillimetre galaxies. 

As a framework for interpreting these data, and for predicting the results of forth- 
coming submm surveys of Swift-derived host samples, we model the expected flux 
and redshift distributions based on luminosity functions of both submm galaxies and 
GRBs, assuming a direct proportionality between the GRB rate density and the global 
star formation rate density. We derive the effects of possible sources of uncertainty in 
these assumptions, including (1) introducing an anticorrelation between GRB rate and 
the global average metallicity, and (2) varying the dust temperature. 

Key words: submillimetre - infrared: galaxies - dust, extinction - galaxies: evolution 
- cosmology: observations - gamma rays: bursts 



1 INTRODUCTION 

There is now strong evidence linking long-duration Gamma 
Ray Bursts (GRBs) with the core-collapse of massive stars 
(e.g. Hjorth et al. 2003a) — and hence, given the short main- 
sequence lifetimes of such stars, with star formation activity. 
Indeed, as tracers of star formation, GRBs hold a number 
of advantages over traditional methods. The high luminos- 
ity of their prompt emission and afterglows enable them to 
be detected, in principle, out to redshifts >10 (in practice 
currently out to z > 6, e.g. Haislip et al. 2006). Their high 
energy emission can furthermore pass unaffected through 
intervening gas and dust — the very conditions one would 
expect to be associated with massive star formation. As the 
outcome of a single stellar event, the luminosity of the GRB 
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ought to be independent of that of its host galaxy, enabling 
localisation of galaxies too faint, dusty or distant to be de- 
tected by traditional means, thus sidestepping many of the 
biases that afflict optical and submm surveys. Furthermore, 
spectroscopy of the bright optical afterglows enables one to 
measure the redshift and other properties of the host galaxy, 
even when direct detection of the galaxy may be infeasible 
(e.g. Berger et al. 2002). 

Once the star-forming properties of a carefully-selected 
subsample of GRB hosts have been established, it should 
be possible to derive the the star formation history of the 
Universe, by measuring the redshift distribution of GRBs. 
A purely GRB-selected galaxy sample should, furthermore, 
represent an unbiased census of the galaxy types reflecting 
their relative contribution to the bulk star formation rate. 
Hitherto it has been difficult to assess the biases affiicting 
the assembly of such samples, factors modulating the GRB 
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rate as a function of redshift (for example, redshift depen- 
dence of density of surrounding medium; metallicity; stellar 
initial mass function (IMF) ; the distribution of jet opening 
angles). Follow-up of samples of bursts detected with Swift 
(Gehrels et al. 2004) shows promise in being able to char- 
acterise and overcome such biases. The BAT (Burst Alert 
Telescope) detector is more sensitive to high-redshift bursts 
than previous missions (e.g. Band 2006), and, due to the 
rapid localisation of bursts via the onboard XRT (X-Ray 
Telescope), ground-based follow up of afterglows (yielding 
information constraining the physical properties of the af- 
terglow, for example redshift, spectral slopes, light curves 
and jet-break times) is much more systematic. 

However, the true proportionality between the global 
GRB and star formation rates has yet to be definitively 
established. It is plausible, for instance, that an otherwise 
direct relation is complicated by dependence on conditions 
local to the GRB. For example, it is thought that metallicity 
could play a role in the GRB formation process (e.g. Fynbo 
et al. 2003; Fruchter et al. 2006). Before one can begin to 
exploit the new GRB catalogue produced by a mission such 
as Swift, it is of the utmost importance to characterise such 
effects. 

One of the most significant contributors to the star for- 
mation rate density at high redshift is the submillimetre 
galaxy (SMG) population (e.g. Small, Ivison & Blain 1997; 
Hughes et al. 1998; Scott et al. 2002; Borys et al. 2003; 
Mortier et al. 2005). Surveys with submm/mm bolome- 
ter arrays such as SCUBA (Submillimetre Common-User 
Bolometer Array) on the James Clerk Maxwell Telescope 
(JCMT), and MAMBO (Max Planck Millimetre Bolometer) 
on the IRAM (Institut de Radioastronomie Millimetrique) 
30m telescope, have revealed a population of galaxies form- 
ing the bulk of their stars in regions optically thick with 
dust, such that most of their extremely high luminosity 
(L ^ 1O^^~^^L0) is emitted in the rest-frame far infrared 
region. Fits to the submm source counts and the FIR back- 
ground imply that the integrated power reprocessed by dust 
contributes a major fraction of the total luminous energy 
emitted throughout cosmic history (Blain et al. 1999). 

Taken together, these facts suggest, all other things be- 
ing equal^, that highly star-forming SMGs ought to yield a 
high rate of GRBs; conversely, that a high fraction of GRB 
host galaxies should turn out to be luminous submillime- 
tre sources. Follow up of GRB hosts in the millimetre and 
submillimetre bands therefore provides one of the most im- 
portant calibrators of the role of GRBs as star formation 
indicators. 

For this reason, in recent years, there have been a num- 
ber of targetted submm studies of GRB hosts, predomi- 
nantly carried out at 850/im using SCUBA on the JCMT 
(Berger et al. 2003, Smith et al. 1999, 2001, Barnard et 
al. 2003). An overview and an analysis is provided by Tan- 
vir et al. (2004). They compared observations made with 
JCMT/SCUBA with model predictions made assuming a di- 
rect proportionality between the GRB rate and the star for- 
mation rate. They discovered that, relative to these predic- 



tions, observations show a deficit of bright (>4mJy) sources. 
Although statistically only marginally significant, this effect 
could, if confirmed, have important ramifications for deriva- 
tions of the global star formation history based on GRB 
surveys. In aU, only three GRB hosts (GRB 000210, ^=0.85; 
GRB 000418, z=l.l2; GRB 010222, ^=1.48) out of 23 with 
850^m RMS values <1.4mJy (Tanvir et al. 2004) have been 
securely detected in the submm. They all lie around the 
3mJy level at 850/im, which, though faint, nevertheless im- 
plies dust-rich (Md ^ lO^^^^Mo), massively star-forming 
(SFR~lOO-lOOOM0yr"^). On the other hand, optical/NIR 
photometry of these hosts reveals stellar populations that 
suffer little dust extinction and have low star formation 
rates (-l-lOMgyr"^) (Gorosabel et al. 2003). This illus- 
trates the importance of submm observations in understand- 
ing the properties of GRB hosts: relying on optical data 
alone, their true significance could easily be overlooked. It 
is thus important to study wider samples of hosts, and es- 
pecially to hunt for additional detections, to test whether 
the existing submm-bright sample is typical or unusual. We 
are in the process of obtaining X-ray observations of the 
submm-detected sample which should settle this question.) 

The present paper pursues these ideas further. The pa- 
per falls into two parts. First, we present new mm and 
submm observations. The former were obtained using the 
MAMB02 bolometer array on the 30m IRAM Pico Veleta 
telescope, as the preliminary part of the first millimetric sur- 
vey explicitly targetting GRB hosts^. The submm data (850 
and 450/im ) were obtained using SCUBA on the JCMT. 
The targets include two of the most extreme GRB hosts 
known (the reddest and the faintest), deliberately selected 
as the most promising candidate submm galaxies. Of par- 
ticular interest, we have obtained deep photometry, at all 
three wavelengths (450/850/1200^m), of the reddest after- 
glow/host found to date, GRB 030115. 

In the second part of this paper, we develop models to 
constrain the relation between GRBs and their host galaxies. 
Using models of the luminosity distribution and evolution of 
submm galaxies, we derive fits to the luminosity function of 
GRBs under the assumption that the GRB rate is a func- 
tion of the star formation rate/far-infrared luminosity of the 
galaxy. Based upon these fits, we estimate the flux distribu- 
tions expected at mm and submm wavelengths, which will 
facilitate comparison between models of the cosmic star for- 
mation history, with future mm and submm surveys of the 
host galaxies of GRBs detected by Swift. 



2 OBSERVATIONS 

2.1 SCUBA/JCMT observations 

The question arises, from previous work, as to whether any 
GRB host galaxies are similar to SMGs For example, the 
three submm-detected hosts all have bluer colours than typi- 
cal of submm galaxies. To address this question, we first used 
JCMT/SCUBA to target a specific host whose properties 
indicate it to be a promising candidate dust-rich, submm 



^ One way in which they may not be equal would be, for exam- 
ple, if a higher than anticipated fraction of SCUBA galaxies were 
powered by AGN 



^ Note that the host of GRB 010222 was scrcndipitously detected 
at 1.2mm during a mm search for its afterglow 



galaxy. GRB 030115 has the reddest optical colours mea- 
sured for a GRB host, implying the presence of a large mass 
of dust, and a (photometric) redshift {z—2.5) placing it near 
the peak of the redshift distribution measured for submm 
galaxies (Chapman et al. 2003). In these respects, this object 
constrasts markedly with the three submm-detected hosts, 
all of which have R~ K < 3 and z < 1.5. 

We obtained new 850pim and 450/im observations of 
GRB 030115 with JCMT/SCUBA on 2005 January 27 and 
28. SCUBA was used in standard photometry mode, with a 
chop throw of 60 arcsec in azimuth. The zenith opacity was 
measured via skydips and the JCMT water vapour monitor, 
and remained within the range 0.065 < r225GHz < 0.08 on 
20050127 and 0.055 < r225GHz < 0.06 on 20050128. Flux 
calibration was obtained from the planets Uranus and Mars 
and several secondary calibrators. Data were reduced both 
manually using the SURF package, and via the ORAC-DR 
pipeline. Additional sky removal was achieved by using off- 
source bolometers to estimate the background. 

GRB 030115 had previously been observed by 
JCMT/SCUBA in Target of Opportunity mode commenc- 
ing 2003 January 18, 3.3 days after the burst, for a total 
of two hours: an upper limit of 6mJy (3ct), at 850^m, was 
reported by Hoge et al. (2003). Since no afterglow was de- 
tected, we can use this measurement as an additional upper 
limit on the submm flux of the host galaxy. We re-reduced 
the archived data in the same manner as described above, 
to find fluxes as reported in Table 1. 



2.2 IRAM-30m/MAMBO2 1.2mm data 

Using the 117-element Max Planck Millimetre Bolometer 
(MAMBO) detector on the Institut de Radiostronomie Mil- 
limetrique (IRAM) 30m Pico Veleta telescope, we obtained 
observations of five GRB hosts between December 2004 and 
April 2005, via pooled (service) observing mode. Selection of 
the targets was designed to improve the redshift distribution 
of the overall submm/mm GRB host sample, in particular 
to try to eliminate a possible bias toward low redshift (see 
Section 3.2.2). All the observed targets lie at 2; > 1, and 
their mean redshift is 2.1. 

Observations were carried out at a wavelength of 1.2mm 
using MAMB02 in On-Off mode. Sky opacity was moni- 
tored frequently by performing skydips; regular pointing and 
focus checks were carried out; and flux calibration was ob- 
tained from standard sources. The data were reduced using 
the NIC software package, which forms part of the GILDAS 
distribution'*. The principles are similar to the SCUBA data 
reduction described above, for example the use of off-source 
bolometers to facilitate sky subtraction. Details of the ob- 
servations, and final fiuxes of the GRB hosts, are reported 
in Table 2. 



2.3 Results 

None of the hosts is detected, either at 1.2mm with MAMBO 
or at 850/450/um with SCUBA. Moreover, the stacked, 
inverse-variance-weighted mean flux of our sample of flve 
new 1.2mm observations is — 0.11±0.27, consistent with a 
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Figure 1. Broad-band SED of the host galaxy of GRB030115, 
showing optical and near-infrared photometry (squares: Lcvan 
et al. 2006) together with SCUBA submm (filled circles) and 
MAMBO mm (unfilled circle) upper limits {2cr). For comparison, 
model SEDs of the template galaxies, the prototype ultralumi- 
nous infrared galaxy Arp220 and the prototype extremely red ob- 
ject HRIO, redshifted to the rest-frame of GRB030115 (assuming 
the photometric redshift z=2.50 derived by Levan et al. (2006)) 
has been plotted. The models were produced by the stellar spec- 
tral evolution code GRASIL of Silva et al. (1998), normalised 
to the NIR/optical flux of 030115. Our submm data indicate 
that GRB 030115 is marginally inconsistent with this Arp220- 
like SED, although an HRlO-like SED cannot be ruled out. Also 
plotted are two isothermal SEDs normalised so as to have the far 
infrared luminosities implied by the optical data and the inferred 
extinction. The dashed line has T=37K, /3=1.5, canonical values 
commonly assumed for submillimctrc galaxies, but is ruled out 
by our submm limits. The dotted line shows that to be consistent 
with the data, hotter dust is required (in this case r=50K). Ob- 
servations in the mid-infrared, for example with Spitzer, would 
be required to fully constrain any hot dust component. 



zero flux for this sample. For comparison, the weighted mean 
850/im flux of the sample discussed in Tanvir et al. (2004) is 
0.93±0.18, which could be interpreted as a true measure of 
the flux of the "typical" GRB host (though as discussed by 
Tanvir et al., the weighted mean carries an "observer bias" in 
that sources with higher fluxes tend to be observed to greater 
depth to attempt to secure detections. The unweighted mean 
of their sample is 0.58±0.36mJy). However it is difficult, 
with such a small sample, to draw any flrm conclusions, and 
we emphasise that investigation of the millimetric properties 
of GRB hosts is ongoing, the present dataset representing 
merely a pilot study. 

Two notable hosts we now discuss individually. 



2.3.1 GRB 030115 

Coadding all the data obtained for GRB 030115, we find flux 
densities 0.0±0.8mJy at 850/im, 7±llmJy at 450/im and 
0.0±0.8mJy at 1.2mm. We note that this GRB was also ob- 
served using MAMBO by Bertoldi et al. (2003), on 20030116 
and 20030118 (i.e. shortly after the burst), to attempt to de- 
tect the afterglow. Their non-detection (0.4±0.9mJy) can, 
again, be combined with our new data, to yield a total 
flux 0.2±0.6mJy. Thus although this host was not detected. 
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Table 1. Summary of JCMT/SCUBA observations of GRB 030115. Zenith opacities are shown at 225GHz, as a range appropriate to 
the time of observations. 

UT date Observation time ''"225GHz 850^m flux 450/im flux 

(min) (uiJy) 



20050127 


40 


0.068-0.070 


3.2±1.9 


71±34 


20050128 


120 


0.055-0.058 


-l.Oitl.l 


-1±12 


20030118 


115 


0.083-0.086 


0.0±1.7 


9±65 



Table 2. Details of our new MAMB02 1.2mm observations of GRB hosts. Zenith opacities are shown as a range. Quoted fluxes are the 
final values obtained by coadding all the datasets 



GRB 



UT Date Observing Time Opacity (r) Final 250GHz flux density 
(min) (mJy) 



GRB 020124 


3.20 


20050201 


60 


0.17-0.19 


0.28±0.60 






20050224 


30 


0.17-0.18 








20050226 


30 


0.25-0.34 








20050327 


30 


0.30-0.35 




GRB 021211 


1.01 


20050223 


30 


0.09-0.10 


0.07±0.53 






20050224 


30 


0.18-0.20 








20050226 


30 


0.25-0.34 








20050327 


30 


0.30-0.35 




GRB 030115 


2.5 


20050223 


25 


0.14-0.15 


0.01±0.76 






20050224 


30 


0.17-0.18 




GRB 030226 


1.98 


20050120 


60 


0.26-0.28 


-0.29±0.66 






20050224 


30 


0.10-0.12 




GRB 030227 


1.6 


20041217 


95 


0.17-0.20 


-0.54±0.53 






20050226 


30 


0.25-0.34 








20050327 


30 


0.24-0.30 





we nevertheless possess reasonably strong upper limits at 
all three wavelengths. (We stress, however, that the 450/ini 
measurement in particular carries a substantial calibration 
uncertainty. Such a deep short-wave limit is rare, but it 
must be used with caution.) The 850/im RMS, in partic- 
ular, would have been easily sufficient to have detected the 
three submm-bright GRB hosts known to date (e.g. Tanvir 
et al. 2004). 

This object is of particular importance in understand- 
ing the host galaxies of GRBs. As noted, it is the reddest 
GRB host observed to date, with a colour R — K ^ 5 (Lcvan 
et al. 2006: hereafter LOG) that qualifies it as an ERO (Ex- 
tremely Red Object). The afterglow, too, is exceptionally 
red {R — K ~ 6), providing further evidence for intrinsic ex- 
tinction. Although no spectroscopic redshift was measured, 
LOG determine a photometric redshift z — 2.5±0.2. Adopting 
the calibration determined by Meurer, Heckman & Calzetti 
(1999) for local starburst galaxies, the extinction implied by 
the rest-frame ultraviolet slope can be used to estimate the 
far infrared luminosity from the observed optical flux, as- 
suming that the absorbed UV photons are reradiated in the 
FIR/submm. The extinction at I6OOA is estimated this way 
to be 5.3, giving a predicted Lfir ~2.4x1O^^L0. Depending 
on the adopted stellar (IMF) initial mass function and star 
formation history, this luminosity implies a formation rate 
of massive stars in the range ~lOO-5OOM0yr^^. Converting 
the luminosity to a predicted submm flux depends upon the 



assumed dust temperature, but for a "typical" ULIRG SED 
(T=40K, 13=1.5), the predicted 850pm flux is ?s2.5mjy. This 
is inconsistent with our measurement. 

In Figure 1 we plot the optical and near infrared pho- 
tometry for this galaxy, together with the new mm and 
submm upper limits. For comparison, we have also plot- 
ted (shifted to the rest frame of GRB 030115) SEDs of 
the canonical ultraluminous infrared galaxy Arp220, and 
the canonical extremely red galaxy HRIO. If this Arp220 
template is normalised to the optical/NIR points, it is 
marginally inconsistent with the submm limits, although 
this represents an extreme case, and the data cannot rule out 
an HRlO-like SED. The dashed line, meanwhile, shows an 
isothermal SED (r=37K, 13=1.5) having the FIR luminosity 
predicted from the optical/NIR flux and inferred extinction. 
Again this is inconsistent with the submm limits, but can be 
accommodated if the dust temperature is increased (dotted 
line) . 

2.3.2 GRB 020124 

Another noteworthy member of the sample is GRB 020124. 
Spectroscopy of its afterglow revealed a high column-density 
(A'^H^IO^^'^) Dampled Lyman Alpha (DLA) system at 
2:=3.2 (Hjorth et al. 2003b). Nonetheless, the host galaxy 
was not detected in deep optical searches with HST, down 
to i? > 29.5 (Berger et al. 2002). One explanation for this 
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faintness is that the host galaxy is dust-rich, and therefore 
plausibly a mm/submm source. Against this interpretation, 
the small extinction {Av < 0.2) inferred from afterglow red- 
dening implies a low gas-dust (Hjorth et al. 2003b) , although 
it is possible that dust in the vicinity of the GRB is destroyed 
by its intense, beamed radiation (Waxman & Draine, 2000) 
hence the line of sight to the afterglow may not be repre- 
sentative of the host galaxy as a whole. Nevertheless, our 
1.2mm upper limit constrains the possible dust mass of the 
host galaxy to < lO^M© (varying inversely with assumed 
dust temperature), which would tend to disfavour a highly 
extinguished host galaxy. 



3 MODELLING THE GRB HOST GALAXY 
SUBMM FLUX DISTRIBUTION 

In this section we explore the degree to which existing and 
future submm observations might be able to constrain the 
efficacy of GRBs as tracers of the star formation rate. Using 
models of the luminosity and redshift distributions of both 
SMGs and GRBs, we can predict the submm flux distribu- 
tion that would be obtained for GRB-selected galaxies. 



3.1 Method and assumptions 

A preliminary calculation along these lines was first per- 
formed by Ramirez- Ruiz, Trentham & Blain (2002). We 
adopt their assumptions concerning the properties and evo- 
lution of SMGs, as described by the models of Blain et 
al. (1999, 2001). Briefly: the submm luminosity function 
is based on the local 60/im luminosity function (Saun- 
ders et al. 1990), with luminosity evolution described by 
$(!/, 2:) = n{z)(f>{L / g{z)) . All submm galaxies are assumed 
to have an identical isothermal dust SED, whose tempera- 
ture is a parameter that can be determined by insisting that 
counts in the mid-infrared are jointly fit. When an emissiv- 
ity index /3=1.5 is assumed, a best fit temperature r=37K 
results from fits to ISO and SCUBA counts (See Blain et al. 
1999 and Blain 2001 for full details.) 

To model the properties and statistics of GRBs, we first 
of all assume that the gamma-ray spectrum is described by 
a Band (1993) function, using values of the spectral indices 
Q=— 1, l3=—2 and cut-off energy _Bo=200keV (rest- frame) . 
We investigate two common parametrizations of the peak 
luminosity distribution: (1) a log- normal luminosity func- 
tion (specified by a mean luminosity Lo and a width a), 
and (2) a Schechter (1976) function (specified by a char- 
acteristic luminosity L» and an index 7). In each case, 
p{z), the (comoving) GRB rate density as a function of 
redshift is derived from the global star formation rate den- 
sity, pgrb(z) ~ ??GRB X ^*(z) — e.g. as calculated from the 
submm models. Similarly, the GRB rate per galaxy is as- 
sumed to be proportional to the galaxy's far-infrared lumi- 
nosity. Initially, we assume that ?7grb — the "efficiency" of 
GRB production — is a constant. However, in general we 
may consider cases where r; is a function of redshift, or of 
host galaxy properties (see below). The parameters Lo, L*, 
a and 7 are then determined, for each star formation history, 
by fitting the flux distribution of long-duration (tgo > 2s) 
GRBs in the BATSE 4B catalogue (Paciesas et al. 1999). 
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Figure 2. Predicted fraction of GRB hosts brighter than a 
given flux density, plotted for four different wavelengths in the 
submm/mni regime. This model displayed here is assumes a 
log-normal high-energy luminosity distribution for long-duration 
GRBs (parameters estimated by fitting to the peak flux distribu- 
tion of the BATSE-4B catalogue), and is calculated for the sensi- 
tivity of Swift/HAT. However, varying these parameters does not 
have a dramatic effect upon the results, at least compared with 
the other uncertainties involved (for example, the SED assumed 
to describe submm galaxies and GRB hosts). 



For this purpose, we adopt a limiting photon flux sensitiv- 
ity 0.27 ph cm~ s~^, corresponding to the median of the 
BATSE sensitivity distribution determined by Guetta, Pi- 
ran & Waxman (2005). 



3.2 Results 

3.2.1 Flux distribution 

In Figure 2, we show the predicted cumulative fraction 
of GRB hosts above a given flux density, for a range of 
submm/mm wavelengths. In this case, we have assumed 
a detector sensitivity appropriate for SwiftfBAT, based 
upon the intercomparison between BAT and BATSE made 
by Band (2006). These results correspond to a log-normal 
GRB luminosity distribution. However, calculations for a 
Schechter function give very similar results. 

In general, the GRB descriptors appear not to affect 
the relative numbers signiflcantly (so long as they represent 
reasonable fits to the number counts). The adopted submm 
properties have a larger effect. In particular, we have as- 
sumed that GRB host galaxies share a common, isothermal 
SED — an SED moreover identical to that of submm galax- 
ies. In reality it is possible that the mean dust temperature 
of GRB hosts is different from the 37K assumed here, and 
that across the sample a distribution of temperatures is to 
be found. Testing the effect of this on our predictions is 
not trivial, since the dust temperature is a parameter of the 
Blain et al. (1999) submm galaxy evolution models. A de- 
tailed refit of the SCUBA counts is beyond the scope of this 
present work, but for now we can obtain a simple indication 
of the effect of varying the temperature by taking note of the 
correlations between the uncertainties in the model param- 
eters of Blain et al. (1999), and rerunning our calculation 
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Figure 3. Effect of varying ttie assumed temperature on two of 
the wavelengtlis plotted in Figure 2. The bold lines correspond to 
the "best fit" temperature as used in Figure 2, while light lines 
illustrate a plausible range of temperatures (30— 50K, top to bot- 
tom). Recall that the submm upper limits on the host of GRB 
030115 imply a temperature toward the upper part of this range. 
In reality, both GRB hosts and the submm galaxy population in 
genernal are likely to exhibit a distribution of temperatures: this 
calculation assumes one SED for all. A hotter than average tem- 
perature for the GRB hosts might explain the paucity of submm 
detections to date. 
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3.2.2 Predicted redshift distribution 

Figure 4 compares the predicted redshift distribution of all 
GRBs (thick solid line) with observed, spectroscopically- 
derived redshifts (light-shaded histogram). Also shown are: 
the variation of the expected fraction of submm-bright hosts 
with redshift, ^n,(S'85o > SmJy) / ■j^n{total) (thick dashed 
line); and the redshift distribution of the submm-observed 
sample (dark-shaded histogram). Some care must be taken 
when interpreting the observational data, since the observed 
distribution is derived from a rather inhomogeneous input 
sample. Not only does the sample consist of bursts detected 
by a range of missions, but, insisting upon spectroscopic 
follow-up inevitably introduces strong biases — for example, 
toward lower-redshift bursts, toward those that are intrinsi- 
cally brighter, or toward those suffering less dust extinction 
from their hosts. 

From the figure it is clear that most of the existing 
submm-observed GRB sample lies at lower redshift than 
the predicted peak in the submm-bright fraction — and in- 
deed the observed peak of the redshift distribution of submm 
galaxies (Chapman et al. 2003). This redshift bias is, there- 
fore, another possible explanation of the lack of submm de- 
tections in the existing host sample. Now, however, afterglow 
redshift determination is more systematic, with the accurate 
localisation provided by the XRT and UVOT instruments on 
board Swift, and rapid ground-based follow-up via a suite of 
robotic and semi-robotic telescopes. Submm/mm follow-up 
of samples resulting from such campaigns is likely to place 
much more secure constraints on the star-forming properties 
of GRB host galaxies than has been possible hitherto. 



using the new parameters. Results for a plausible tempera- 
ture range are shown in Figure 3. Although the uncertainty 
is probably somewhat exaggerated (because it assumes all 
submm galaxies are affected in the same way) the effect is 
much larger than that due to any uncertainty in the GRB 
properties. The submm SEDs of GRB hosts are ill-enough 
constrained that it seems plausible that a hotter than aver- 
age dust temperature — as our limits for GRB 030115 would 
imply — could account for the paucity of detections to date. 
At 1.2mm, we would expect to have detected ~10-15 
percent of a sample with RMS ~0.5-0.6mJy, assuming the 
T=37K model. This could increase to as much as --^20-25 
percent if the temperature were permitted to be as low as 
30K, but would become negligible at temperatures as high 
as 50K. The predicted flux distributions enable us to calcu- 
late the average flux density of a large sample, to compare 
with the coadded (detected plus non-detected) fluxes from 
observations. The 37K model predicts < 5*850 >~ O.SmJy, 
< Si. 2 >^ 0.4mJy. Recall that Tanvir et al. (2004) find 
stacked fluxes 0.93±0.18 (weighted mean) and 0.58±0.36 
(unweighted), both consistent with this prediction. Our 
1.2mm mean, on the other hand, is marginally inconsistent 
with the prediction. Our sample of five is, however, too small 
to confirm or reject any of the models, but continued study 
of homogeneously-selected host samples should improve the 
constraints. Ultimately, greater sensitivity will be attained 
by taking advantage of forthcoming facilities such as ALMA 
(or even existing facilities such as Spitzer): then it will be 
possible to reach limits deep enough to discriminate between 
models. 



3.2.3 Dependence of GRB rate on metallicity 

One important factor that might ultimately mitigate against 
the formation of GRBs in dust-rich galaxies is the role that 
metallicity is thought to play (e.g. Fruchter et al. 2006). Ac- 
cording to the "coUapsar" model (MacFadyen & Woosley, 
1999), a low metal abundance allows the progenitor to retain 
a high mass and angular momentum, favouring the produc- 
tion of a black hole and accretion disk. High detection rates 
of Lyman-a emission from GRB hosts (Fynbo et al. 2003) 
could be taken as evidence that these systems are indeed 
metal-poor. If this metallicity dependence is correct, it hold 
consequences both negative and positive for the use of GRBs 
as star formation indicators. Whilst, on the one hand, com- 
plicating the conversion between GRB and star formation 
rate, it suggests that GRBs may instead be the ideal means 
of pinpointing metal-poor galaxies — in particular low-mass, 
unenriched systems at the highest redshifts which are most 
likely to be missed in other surveys. 

It is therefore important to explore the possible effects 
of a metallicity dependence. To do so, we place a redshift 
dependence on the GRB rate density — SFR density conver- 
sion factor ?7grb(z). The evolution of the average metallicity 
with redshift is given by the submm galaxy models (Blain et 
al. 1999). This is converted to a relative efficiency of GRB 
production using an ad hoc recipe — which ultimately may 
be unrealistic, but, in the absence of any compelling obser- 
vational or theoretical guidelines, it serves amply to illus- 
trate the effects. A sample predicted redshift distribution 
is plotted in Figure 4 (thin curve). As expected, the peak 
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Figure 4. Predicted (curves) and observed (histograms) redshift 
distribution of GRBs. The thick, continuous curve is the predic- 
tion from the basic submm galaxy model, assuming a detector 
sensitivity appropriate for Swift/H AT. The thick dashed line il- 
lustrates the variation with redshift of the fraction of submm- 
bright {Ssso > 3mjy) hosts. The light-shaded histogram shows 
the distribution of all GRBs with spectroscopically-confirmed red- 
shifts, while the dark-shaded histogram shows the subset of this 
sample for which sensitive submm photometry has been carried 
out. Finally, the light continuous curve is a variation on the basic 
model incorporating a dependence of the GRB formation effi- 
ciency on the global average metal abundance. The peak is, as 
one would expect, shifted to higher redshift, resulting in a dis- 
tribution more consistent with that emerging from spectroscopic 
follow-up of Swift GRBs. (N.B. the vertical scale is appropriate 
for the dashed line (submm- bright fraction): the remaining curves 
and histograms are scaled arbitrarily.) 



is shifted toward higher z where the average abundance of 
heavy elements is smaller. However, without separately en- 
coding galaxy-to-galaxy variations in the metal abundance, 
the efTects on the submm flux distribution are small. 

This calculation is, we emphasise, only illustrative at 
present. For example, it may be more appropriate, for GRB 
hosts, to consider metallicities traced by optical galaxy sur- 
veys, rather than submm surveys as used here. As the red- 
shift distribution of GRBs becomes more fully sampled, (for 
example via spectroscopic follow-up of large samples of Swift 
bursts), it will soon be possible to place constraints on a 
wider range of models in this way. 



4 SUMMARY 

Following from the previous study of Tanvir et al. (2004), 
we have further investigated the millimetre/submillimetre 
properties of the host galaxies of GRBs, in order to char- 
acterise the efficacy of GRBs as star formation indicators. 
Specific increments over the T04 study include: (1) we have 
conducted the first survey of GRB hosts at millimetric wave- 
lengths, with the MAMB02 bolometer array on the IRAM 
30m Pico Veleta telescope. None of these targets was de- 
tected, down to an average RMS ~0.6mJy at 1.2mm; (2) we 
obtained deep submm photometry of GRB030115, whose 
high intrinsic extinction inferred from its optical/NIR spec- 
tral slope make it a promising candidate submm galaxy. 



Despite its ERO-like optical colours, however, this galaxy 
is not detected in the mm/submm, to deep limits at 850/im 
(cr=0.8mjy) and 450/im (a^llmjy); (3) we have modelled 
the redshift and flux distribution of GRB hosts, assum- 
ing a link between GRBs and the submm galaxy popula- 
tion. A novelty of these models is that they take account 
of the metallicity bias widely proposed to affect the GRB- 
to-star formation rate conversion. As such they potentially 
have much wider applicability than the derivation of submm 
properties, and we will further develop these ideas in future 
publications (Priddey et al., in prep.). 

The non-detection of GRB 030115 is revealing. One 
might contrast this result with the three GRBs that do 
possess submm detections, for their optical/NIR colours are 
much bluer. The broadband spectrum of the GRB 030115 
host is inconsistent with the SED of an extremely luminous 
infrared galaxy such as Arp220 or with a cool, isothermal 
model, but hotter dust (>50K), or template SEDs of other 
submm-luminous galaxies, cannot be ruled out. Observation 
in the mid-infrared with missions such as Spitzer should also 
be able to constrain any hot dust component too faint to be 
seen in the submm. 

We emphasise that this work is ongoing: in the immi- 
nent future we will be able to draw upon larger, post-Swift 
samples of GRBs to ensure a uniform sample selection — 
enabling, for example, a more uniform redshift distribution. 
For the moment, it seems that the trend of a low submm 
detection rate of GRB hosts, seen in previous surveys, is 
maintained. 

What are the implications of a low mm/submm detec- 
tion of GRB hosts? We have shown that there is sufficient 
uncertainty in models and underlying assumptions, as yet 
poorly constrained by observation (for example the adopted 
dust temperature) that a correlation between massive, dust- 
enshrouded star formation and GRB production cannot be 
firmly ruled out. Sample selection biases (e.g. against high 
redshift and highly extinguished bursts) are also likely to 
have played a significant role in previous studies. Our mod- 
els indicate that redshift bias in particular could account 
for the lack of detections within existing surveys. The new 
observations reported here (5 hosts all at z > 1, one highly 
extinguished host and one extremely faint host) were taken 
in part to alleviate such problems. Prior to ALMA, observa- 
tions of consistently foUowed-up samples with existing facil- 
ities (e.g. IRAM-30m, APEX) must be made to enable us to 
make further progress in exploring these effects. The capa- 
bilities of Swift, combined with efficient ground-based follow 
up, show promise in being able to yield such a sample. 



5 ACKNOWLEDGMENTS 

RP thanks PPARC for support. We are grateful to those ob- 
servers and support staff at JCMT and IRAM Pico Veleta 
who carried out observations for us in service mode, and 
we thank the anonymous referee for constructive comments. 
The JCMT is operated by the Joint Astronomy Centre in 
Hilo, Hawai'i, on behalf of the parent organisations: the Par- 
ticle Physics and Astronomy Research Council in the United 
Kingdom, the National Research Council of Canada and the 
Netherlands Organisation for Scientific Research. 



8 Priddey et al. 

REFERENCES 

Band D.L., et al., 1993, ApJ, 413, 281 

Band D.L., 2006, astro-ph/0602267 

Barnard V.E., et al., 2003, MNRAS, 338, 1 

Berger E., et al., 2002, ApJ, 581, 981 

Berger E., Cowie L.L., Kulkarni S.R., Frail D.A., Aussel 
H., Barger A. J., 2003, ApJ, 588, 99 

Bertoldi F., Frail D.A., Berger E., Menten K.M., Kulkarni 
S., 2003, GCN circ. 1835 

Blain A.W., Small I., Ivison R.J., Kneib J.-R, 1999, MN- 
RAS, 302, 632 

Blain A.W., 2001, in Starburst Galaxies: Near and Far, 
Proceedings of a Workshop Held at Ringberg Castle, Tac- 
coni L. & Lutz D. (Eds.), Heidelberg: Springer- Verlag, 
p.303 

Borys C, Chapman S., Halpern M., Scott D., 2003, MN- 
RAS, 344, 385 

Chapman S.C., Blain A.W., Ivison R.J., Small I.R., 2003, 
Nature, 422, 695 

Frail D.A., et al., 2002, ApJ, 565, 829 

Fruchter A.S. et al., 2006, Nature, submitted 

Fynbo J.P.U., et al., 2003, A&A, 406, L63 

Hjorth J., et al., 2003a, Nature, 423, 847 

Hjorth J., et al., 2003b, ApJ, 597, 699 

Hoge J.C., Stevens J. A., Moriarty-Schieven G., Tilanus 
R.P.J. , 2003, GCN circ. 1832 

Hughes D., et al., 1998, Nature, 394, 241 

Gorosabel J., et al., A&A, 409, 123 

Gehrels N. et al., 2004, ApJ, 611, 1005 

Guetta D., Piran T. & Waxman E., 2005, ApJ, 619, 412 

Haislip J.B., et al., 2006, Nature, 440, 181 

Levan A. et al., 2006, ApJ, in press (ApJ preprint 
doi:10.1086/'503595') 

MacFadyen A.I. & Woosley S.E., 1999, ApJ, 524, 262 

Meurer G.R., Heckman T.M. & Calzetti D., 1999, ApJ, 
521, 64 

Mortier A.M.J. , et al., 2005, MNRAS, 363, 563 

Paciesas W., et al., 1999, ApJS, 122, 465 

Ramirez-Ruiz E., Trentham N., Blain A.W., 2002, MN- 
RAS, 329, 465 

Saunders W., Rowan-Robinson M., Lawrence A., Efs- 
tathiou G., Kaiser N., Ellis R.S., Frenk C.S., 1990, MN- 
RAS, 242, 318 

Schechter P., 1976, ApJ, 203, 297 

Scott S.E., et aL, 2002, MNRAS, 331, 817 

Silva L., Granato G.L., Bressan A., Danese L., 1998, ApJ, 
509, 103 

Small I., Ivison R.J., Blain A.W., 1997, ApJ, 490, L5 

Smith LA., et al., 1999, A&A, 347, 92 

Smith LA., Tilanus R.P.J. , Wijers R.A.M.J., Tanvir N., 
Vreeswijk P., Rol E., Kouveliotou C, 2001, A&A, 380, 81 

Tanvir N.R., et al., 2004, MNRAS, 352, 1073 

Waxman E. & Draine B.T., 2000, ApJ, 537, 796 



